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Abstract 

Using the combined CLEO II and CLEO II. V data sets of 9.1 fb _1 at the T(4S'), we measure 
properties of ip mesons produced directly from decays of the B meson, where ll B" denotes an 
admixture of B + , B , B°, and B°, and "•0" denotes either J/^{1S) or ip(2S). We report first 
measurements of ip polarization in B — > V(clirect)X: awjm = — 0.30to'og ± 0.04 and a^^s) = 
— 0.45^o'i9 i 0.04. We also report improved measurements of the momentum distributions of 
ip produced directly from B decays, correcting for measurement smearing. Finally, we report 
measurements of the inclusive branching fraction for B — > tpX and B — > x&X. 
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Inclusive production of ip is currently understood in the framework of Non-Relativistic 
QCD (NRQCD) effective field theory In 1995, measurements of prompt i)j produc- 
tion at the Tevatron 0, [| ruled out the then-dominant Color Singlet Model (CSM); in 
contrast, NRQCD calculations j3J could accommodate the relatively large production rate. 
However, measurements of the polarization of these prompt ip || deviated from NRQCD 
calculations at high py. The precision of these calculations is limited by the knowledge 
of the process-independent, long-distance matrix elements (LDME's), which also appear in 
NRQCD calculations of ip production in B decays. The polarization of ip produced from B 
decays [^|, [7]] is sensitive to the color-octet LDME's; however, these calculations have been 
done only to leading order (LO). The momentum distribution of ip produced in B decays 
||. [|, |10| is also sensitive to the dominant color octet terms, particularly at low p^. Addi- 
tionally, the low-momentum region would also be affected by the existence of an intrinsic 



charm component in B mesons [|ll|]. Finally, the inclusive branching fraction B(B — > ifiX) 
H [12], [H| constrains a sum of LDME's. This Letter reports measurements of these three 
properties of if: production in B decays, which could significantly improve the knowledge of 
the non-perturbative parameters of NRQCD. 

Our analysis |14| is based on 9.7 million BB events (9.1fb _1 ) produced on the Y(45) 
resonance. Additionally, 4.4 fb -1 of data collected slightly below the T(45) resonance were 
used to subtract the small (~ 2%) contribution of continuum (e + e~ — > qq, q G {u, d, c, s}) 
ifj(lS) production. The e + e~ collisions were delivered by the Cornell Electron Stoarge Ring 
(CESR) and detected with two configurations of the CLEO detector, CLEO II |L5| and 
CLEO II. V 0. 

We select events that have spherical energy distributions and are likely to be hadronic. 
We reconstruct ijj candidates in the dilepton modes ijj — > and -0 — > e + (7)^(7). 

The selection criteria were chosen with a goal of high detection efficiency. In the di-muon 
channel, at least one of the muon candidates must penetrate at least 3 interaction lengths 
into the iron of the solenoid return yoke; if only one candidate satisfies this, then the other 
candidate must leave a shower in the crystal calorimeter which is consistent with that of a 
minimum ionizing particle. In the di-electron channel, we use shower information from the 
crystal calorimeter and measurements of specific ionization from the drift chamber to identify 
electron candidates. We also attempt to recover up to one Bremsstrahlung photon for each 
electron candidate. To do this, we select the most collinear shower within a five-degree cone 
around the initial electron direction; furthermore, the shower must not be associated with 
any track, and, when combined with any other shower in the event, must not result in an 
invariant mass consistent with a 7r°. The PDG 2001 [17j branching fractions are used to 
combine results from the electron and muon channels, except for £>(-0(25) — > which 
we assume by lepton universality to be equal to £>(-0(25) — > e + e~), with an uncertainty of 
Vo of itself; this is consistent with recent measurements HI! 



About 30% of ip(lS) from B decays have intermediate parents of ij)(2S) or Xci [17 



Our measurements of directly produced i/j(1S) are obtained by subtracting the contribu- 
tions of these "feed-down" i(j(1S) from the inclusive '0(15') sample. For every event with 
a ip(lS) candidate within MeV of the nominal mass, we search for these intermediate 
parents through the decay chains Xci —* ^ , (1<S')7 an d '0(25') — > -0(15)7t + 7t~. We reconstruct 
(Mg+g- n + n - — Mi+i-) and (M^+£- 7 — Mi+£~), which have better resolution than the recon- 
structed "0(25) and Xci invariant masses themselves. In the ip(2S) —>■ if)(lS)iT + 7i~ decay 
chain, we reduce low-momentum pion background by requiring M n + n - > 0.45 GeV, which 
has a efficiency of about 85%, from Monte Carlo simulation. This decay mode is also used 
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to improve the statistics in our measurements of the inclusive branching fraction and the 
ip(2S) momentum distribution in B — > r ip(2S)X. We do not reconstruct the related decay 
ip(2S) — > ifj(lS)7r Ti , and argue that the properties of ij)(lS) from this decay are identical 
to those of ip(lS) from ip(2S) — > t(j(lS)n + n~; the kinematic difference in the momentum 
distribution is small compared to the experimental resolution, and the isospin state of the 
7T7T state has no bearing on the polarization of the ip(\S). 

The CLEO Monte Carlo simulation, based on GEANT [BJ, is used to obtain the invariant 
mass lineshape for signal events and to estimate the detection efficiency In these simulated 
signal events, one of the B mesons decays via one of the decay chains listed above. For each 
decay chain, we generate two samples of events; one with all ip longitudinally polarized, the 
other with all ip transversely polarized. We find that the detection efficiency varies slightly 
SIS db function of ip momentum and polarization. 

The procedure and results for the inclusive branching fraction and momentum distribu- 
tion measurements are as follows. We divide the data into partitions in p^, the momentum 
of the t/> candidate, using a binsize of 0.1 GeV/c. For each partition, the invariant mass 
distribution of ip candidates is fit to a sum of a signal lineshape, obtained from the Monte 
Carlo simulation, and a cubic polynomial background. The average x 2 of the fits is consistent 
with the number of degrees of freedom, thus justifying our choice of the above parametriza- 
tion. We repeat this procedure using signal Monte Carlo events, binning in generated ip 
momentum, to obtain detection efficiencies as a function of The data is then corrected 
for detection efficiency bin by bin; this minimizes the effect of any discrepancy between the 
true p^p distribution and that generated by the Monte Carlo simulation. Similarly, we fit 
the invariant mass distributions of {M^+i- 1T + ll - — M^+i-) and (M^+£- 7 — Mg+i-), to extract 
efficiency-corrected yields of feed-down ip(lS). We thus obtain momentum distributions of 
i/j(1S) and i/)(2S) which have been corrected for detection efficiency, tp(lS) feed-down, and 
continuum background. The yields are then normalized by n^x B(if) —>■ £ + £~(n + n~)), where 
ub is the number of B and B mesons in the data; the uncertainties in these quantities are 
reflected in our results as an overall scale factor error. Inclusive branching fractions are 
obtained by summing the normalized momentum distributions over all bins. Finally, the 
Monte Carlo simulation is used to obtain a matrix which correlates the momentum of the 
xjj as generated to the momentum as measured; by inverting this matrix and applying it 
to the observed momentum distribution, we are able to deconvolve the effects of detector 
measurement smearing from the distribution. 

We investigate the possible sources of systematic error; for each source, we make an 
appropriate modification to the measurement procedure and observe the deviation of the 
resulting yield relative to the nominal procedure. The deviations are then combined to 
obtain final systematic errors. The sources of error are grouped as follows: (1) Monte Carlo 
simulation of track and shower finding, electron and muon identification, if; polarization, 
global event and kinematic cuts, (2) invariant mass fit procedure; (3) branching fractions of 
unmeasured modes; and (4) overall scale factor. 

The results for the inclusive branching fractions are given in Table | and the momentum 
distributions are shown in Figure Our branching fraction results are consistent with pre- 
vious published measurements [f20| as well as preliminary measurements [21J , and are limited 
by systematic errors. The combined error for B(B — > if; (IS) (direct) X) is smaller than the 
error of the PDG 2001 average [117] by a factor of two. However, the theoretical uncertainties 
in the NRQCD calculations of the branching fractions are such that the improved accuracy 
of this measurement is unlikely to further constrain the NRQCD LDME's. The momentum 
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distributions reported here are the first to subtract the measured distributions of feed-down 
and continuum ip, correct for detector measurement smearing, and analyze systematic errors 
for each bin individually. Figure 1(b) is also the first to show the momentum distribution of 
multibody (> 3-body) decays in B — > tfj(2S)X production; these decays account for much 
of the total B — > ip{2S)X production, as is also the case with B — > i/j(1S)X [[J. It should 
be possible to update previous phenomenological studies of the if) momentum distribution 
with these improved measurements. 



Decay 




Branching Fraction (%) 


B - 


+ if>(lS)X 




1.121 ± 0.013 ± 0.040 ± 0.013 


B - 


-> ^(15) (direct) A" 


0.813 ±0.017 ±0.036 ±0.010 


B - 


-> XdX -» 


^{1S)X 


0.119 ±0.008 ±0.009 ±0.001 


B - 


^XdX 




0.435 ± 0.029 ± 0.031 ± 0.026 


B - 


-> ^(2S)X 


-> ip(lS)X 0.189 ± 0.010 ± 0.018 ± 0.002 


B - 


-» ip(2S)X 




0.316 ± 0.014 ± 0.023 ± 0.016 



TABLE I: Inclusive branching fraction results. The errors shown are (in order) statistical, sys- 
tematic, and due to an overall scale factor uncertainty. 

The polarization parameter a is equal to (±1, 0, —1) for a population of (transversely, 
randomly, longitudinally) polarized if). For ifi — > £ + £~ decays, it is determined experimentally 
by measuring the decay angle 9, which is defined as the angle between the £ + direction in 
the ip rest frame and the ip direction in the B rest frame. The cos 9 distribution for a 
population of ijj is proportional to (1 + a cos 2 9). The angular distribution is obtained in a 
similar manner as the momentum distributions: the dataset is partitioned into 5 equal bins 
in cos (9 between —1 and 1; for each partition, we fit the invariant mass distribution to find 
the signal yield. In addition to measuring the polarization of direct ift(lS) and tp(2S) for all 
momenta, we also extract a^is) m 3 coarse momentum bins. 

At CLEO, B mesons are produced with a small boost in the T(AS) (lab) frame, the 
direction of which is unknown. The boost of the B results in a smeared measurement of 
cos 9; directly fitting for a using the measured cos 6* distribution would yield a biased result. 
However, this kinematic smearing is accurately modeled by the Monte Carlo simulation. Our 
procedure is to generate Monte Carlo events in two sets; one with all x/j generated longitudi- 
nally, the other with all transverse. The measured cos 9 distribution from the data is then fit 
to a sum of the reconstructed cos 9 distributions from the polarized Monte Carlo sets. This 
procedure correctly accounts for both the boost smearing and detection efficiency. Since the 
efficiency also depends on p^, we must ensure that the Monte Carlo distributions of generated 
p^p match those of Figure |l|; this is accomplished through a rejection technique. Because the 
observed cos 9 distributions are not directly corrected for detection efficiency, the observed 
feed-down distributions are corrected only for the efficiency of detecting the additional parti- 
cles needed to reconstruct the i()(2S) or Xci- The final feed-down and continuum-subtracted 
angular distributions are shown in Figure The systematic error study included the previ- 
ously mentioned sources of bias; additionally, we have investigated the possible systematic 
error arising for the methods for feed-down subtraction and fitting for a. 

The final polarization results are listed in Table f|; these are the first results for the 
polarization of ifi(lS) and ip(2S) from B — > -0(direct)X. For comparison, we measure 
a = —0.35 ± 0.03 (statistical error only) for ip(lS) from B — > if) (IS) (all) X . Our result 



5 



1.6 
1.4 
1.2 
1.0 



o*0.8 



2030402-003 

i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i 

(a)V(1S) 
• Data Systematic Error Range 



8-^^^(1270) 

B-^^K*(892) 
Sum of Above 



B—^tpJt 




0.4 0.8 1.2 1.6 2.0 
V Momentum (GeV/c) 



FIG. 1: Momentum distributions of (a) ipO-^) an d (b) ^(250 produced directly from B decays. 
There is an additional overall scale uncertainty of 1.2% for ip(lS) and 5.1% for ip(2S) which is not 
depicted in the plots. The histograms show the contributions of two-body B — > tpX decays, where 
the lineshapes are obtained from Monte Carlo simulation and the normalizations are from previous 
determinations of exclusive branching fractions [17, p2|, 23]. 



for a^(is) is about 4 standard deviations from zero; this measurement therefore strongly 
disfavors the color evaporation model of charmonium production P% ], which predicts zero 
net polarization, independent of the production mechanism. When next-to-leading-order 
calculations become available, these results also have the potential to significantly constrain 
the long-distance matrix elements of NRQCD. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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FIG. 2: Decay angle distributions of (a) ip(lS) and (b) ip(2S) from B — > ^(direct)X, summed over 
all pw,. The points represent the data, showing la statistical errors. In both figures, the fit result 
(solid histogram) is the sum of a longitudinal component (dashed histogram) and a transverse 
component. 
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TABLE II: Polarization of ip(lS) and ip(2S) from i? — > ^(direct)X over the full momentum 
range (top two values) and for ip(lS) in three momentum ranges. The errors are statistical and 
systematic. 
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